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e-MERLIN resolves Betelgeuse at A 5 cm: hotspots at 5 R* 
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ABSTRACT 

Convection, pulsation and magnetic fields have all been suggested as mechanisms for the 
transport of mass and energy from the optical photosphere of red supergiants, out to the region 
where the stellar wind is launched. We imaged the red supergiant Betelgeuse at 0.06-0.18 arc- 
sec resolution, using e-MERLIN at 5.5-6.0 GHz, with a sensitivity of ~10 fiJy beam -1 . Most 
of the radio emission comes from within an ellipse (0.235x0.218) arcsec 2 (~5x the optical 
radius), with a flux density of 1.62 mJy, giving an average brightness temperature ~1250 K. 
This radio photosphere contains two hotspots of 0.53 and 0.79 mJy beam -1 , separated by 90 
milli-arcsec, with brightness temperatures 5400 + 600 K and 3800 + 500 K. Similar hotspots, 
at more than double the distance from the photosphere of those seen in any other regime, were 
detected by the less-sensitive 'old' MERLIN in 1992, 1995 and 1996 and many exceed the 
photospheric temperature of 3600 K. Such brightness temperatures are high enough to em- 
anate from pockets of chromospheric plasma. Other possibilities include local shock heating, 
the convective dredge-up of hot material or exceptionally cool, low density regions, transpar- 
ent down to the hottest layer at ~40 milliarcsec radius. We also detect an arc 0.2-0.3 arcsec 
to the SW, brightness temperature ~ 150 K, in a similar direction to extensions seen on both 
smaller and larger scales in the infra-red and in CO at mm wavelengths. These preliminary 
results will be followed by further e-MERLIN, VLA and ALMA observations to help resolve 
the problem of mass elevation from 1 to 10 in red supergiants. 
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1 INTRODUCTION 

Betelgeuse (a Ori), type M2I ab, is the closest re d supergiant, with 
an initial mass of 15-20 M s dDolan et alj|2008h . The asymmetric 
SI and S2 CO outflows have average velocities and radius limits 
of 9.8 km s~*, 4 arcsec and 14.3 km s -1 , ~17 arcsec, respectively 
l lO'GonnanetaLll2012l) . The mass-loss rate M of Betelgeuse is 
~10~ 6 Ms yr ~\ estimated bv lLe Bertre et al.l d2012h from a 0.24-pc 
HI shell. Its proper motion is t owards the NE and a bow sh ock was 
identified in the infra-red (IR) ( Noriega -Crespo et al . 1997) and HI 
and mapped using Herschel jPecin et alj|2012h . which shows arcs 
6-7 arcmin NE of the star. W e adopt the 2.2 /an photospheric ra- 
dius -22.5 milli-arcsec (mas) dDvck et alJl996t|Perrin et al.l2004) 



as R+, at a distance of ~ 197 pc dHarper et al 



2008 



Mass loss from cool, evolved stars (e.g. VX Sgr) which form 
copious dust at ~5— 10 R*, is thought to be initiated by pulsations, 
followed by ra diation pressu re on dust accelerating the wind to es- 
cape velocity dBowerjTl988h . However, the requirements for this 
mechanism do not seem to be present in early M-type stars like 
Betelgeuse. It shows weaker optical variability (Percy et al. 1996) 
than the RSG which form more dust. IVerhoelst et alj (12006) found 
signs of metal oxide nucleation at ~1.5 R* but silicates have only 



been f ound at radii >500 ma s Ce.g. lTatebe et al.l2"0 07: Skinn eret al.l 
1 19971) . IPupree et al.l d 19871) measured periodicities in optical and 
ultraviolet fluxes probably tracing pulsatational shock waves which 
could help to initiate mass outflow. Other potential mass-loss mech- 
anisms i nclude localised events related to conv e ction (e.g. as mod- 
elled by IChiavassa et alj|201p|) . lAuriere et all d2010l) measured a 
longitudinal magnetic field of 0.5 -1.5 G, leading to mode ls for 
magnetic promotion of mass loss dThirumalai & Hevlll2012l) . See 
the Proceedings of the Betelgeuse 2012 workshop, eds. Kervella et 
al., in prep., for more comprehensive reviews of mass loss studies. 



VLA observations by iLim et all dl998l) at 43, 22, 15, 8 and 
5 GHz showed the measured radius increasing from 2-7 The 
opacity increases with decreasing frequency, so lower frequency 
images sample layers at larger radii, found to be cooler. Over 
this range of R*, the radio brightness temperature 7b falls from 
~3540 - 1370 K and the photo sphere becomes c ool enough for 
dust to begin to form. However, Uitenbroe k et al.l (1 19981) detected 
UV continuum and lines formed at >5000 K, at radii of up to 2.8 
and 6 R+, respectively. Absoiption lines reveal a rotation speed of 
5 km s -1 and they suggest that a bright spot in the SW defines the 
orientation of the rotational axis, at a position angle of 55°, giving a 
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deprojected rotation period of 25 yr at 14° yr -1 . The chromospheric 
tracer Mgll sho ws an expansion velocity o f 10 km s -1 , detected out 
to almost 9 R+ (Gilliland & Dupree 1996). Asymmetric ab sorption 
lines such as Ha JWevmannlll962l iBagnulo et al. 1120031) suggest 
outflo w of chromosph eric material at 5-7 km s . 

iLim et "aT] dl998h deduced that the chromospheric gas must be 
at least a thousand times less abundant than the material dominating 
radio emission, ruling out global heating as a mass -loss mechanism, 
leavin g convection as a plausible alternative. Ha rper & Brownl 
(2006) studied multi-epoch spatially resolved HST spectra and, by 
measuring the electron density and comparing it to 22-GHz obser- 
vations, confirmed the small filling factor. Chromospheric emission 
lines suggested rotation about an axis at a position angle of 65° at a 
velocity not exceeding ~ 5 km s -1 , suggesting that chromospheric 
material is not in solid-body rotation at R > 75 mas. 

We observed Betelgeuse using e-MERLIN at 5.5-6.0 GHz 
(/l«5.2 cm), providing ~ 4 resolution elements across the stellar di- 
ameter. These are the first well-resolved images at cm wavelengths, 
investigating whether the irregularities seen at 43 GHz (A7 mm) 
persist at larger radii. We describe the observations in Section [2] 
and results in Section[3] discussing these in Sectionf4] We suggest 
future work in Section|5]and conclude with a summary in Section[6] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 2012 e-MERLIN observations 

Betelgeuse was observed in the first semester of e-MERLIN open 
time on 13-15 July 2012, with a bandwidth of 0.512 GHz centred 
on 5.75 GHz. Seven antennas were used including the 75-m Lovell 
telescope providing baselines from 11-217 km (90-3910 kA). The 
data were processed in dual polarization, using 4x128 MHz spec- 
tral windows, each divided into 64 2-MHz channels. The point-like 
QSO 0551+0829, separation l.°5, was used as the phase reference 
on a cycle of 7:3 min. OQ208 was used as the bandpass and flux 
density calibrator. Calibration is summarized in Appendix A. The 
flux scale is accurate to ~10% and the astrometric uncertainty aris- 
ing from phase referencing is 16 mas. The calibrated and edited 
Betelgeuse data comprised 4-8 hr per antenna, spread over 10.5 hr, 
with an average bandwidth of 400 MHz. 

We made two maps after self-calibration, using different 
weighting schemes (see Appendix A). We maximised sensitivity 
using a 3500-k/l taper and a circular 180-mas (FWHM) restoring 
beam. This produced a noise cr ims 0.027 mJy beam -1 . We max- 
imised resolution by applying partial uniform weighting (ROBUST 
0.75) and using a restoring beam of (80 x 60) mas 2 at position angle 
(PA) 143° (close to the natural beam shape), giving cr rms 0.009 mJy 
beam -1 . Despite the lower noise per beam, the minimum brightness 
temperature sensitivity limit is 225% higher with this weighting. 
In order to ensure artifacts were not introduced, we made images 
with completely natural weighting (used for astrometry) and with 
the above weightings, prior to self-calibration. The high-resolution 
dirty map has <r lms of 0.03 mJy beam -1 , with peaks of 0.81 and 0.63 
mJy beam -1 , similar to the features described in Section |3~3| 



2.2 Archival data 

We retrieved MERLIN Betel geuse data obs erved at 4.994 GHz on 
1995-06-24 and 1996-11-03 jMorrisll200lh using 16 MHz band- 
width, and VLA data observed at 4.885 GHz on 1996-10-21 us- 
ing 100 MHz bandwidth. These data sets were reduced using stan- 
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Figure 1. The main panel shows the 2012 e-MERLIN colour image of 
Betelgeuse optimised for sensitivity to extended structure, using a beam 
size of 1 80 mas, cr nIls = 0.027 mJy beam - 1 . The contours of the image opti- 
mised for resolution are overlaid, with the beam size at lower left, the noise 
o- rras = 0.009 mJy beam -1 , and contour levels at (-1, 1, 2, 4, 8, 16) X 0.027 
mJy beam -1 . The small, solid yellow circle is the size of the 45-mas diam- 
eter photosphere, centred on the low-re solution peak, an d the large, dashed 
circle shows the 310-mas, 5-GHz disc iLim et al . 1998). The insert shows 
the SW arc after subtracting the central star, on the same scale. 



dard techniques dPiamond et ai]|2003l lGreisenlll994l) . The MER- 
LIN data were imaged with natural weighting (beam size (85 x 55) 
mas 2 ) and no cleaning, giving <r rms of 0.15 and 0.12 mJy beam -1 
in 1995 and 1996, respectively. The 1996 MERLIN and VLA data 
were combined to give a resolution of 200 mas and cr rms of 0.058 
mJy beam -1 ; the VLA-only image had 400-mas resolution. 



3 RESULTS 

3.1 Extended structure seen in 2012 

Figure Q] shows the colour image made with optimum sensitivity 
to extended structure, using a 180-mas restoring beam. The peak 
is 1.116 mJy beam -1 . The 3cr rms boundary (0.081 mJy beam -1 
at this resolution), represented by the limit of red shading, has a 
quite irregular outline, with a maximum extent of 550 mas approx- 
imately N-S. The shortest, radial separation between the peak and 
the 30-nns limit is about 190 mas, to the NNE. We fitted a 2D ellipti- 
cal Gaussian component to this image, yielding a total flux density 
1.619 ± 0.057 mJy at position 05:55:10.3274 +07:24:25.514, with 
a noise-based uncertainty of 4 mas. The component had major and 
minor axes (0 ma j, 6 min ) of (235x218) mas 2 , uncertainty tr e = 6 mas, 
PA 1 15 + 9° (significantly different from the beam PA of 143°). 

The brightness temperature of a component of flux density S , 
area A, is given in Kelvin by 
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0.052 m. The area of a Gaussian 
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Figure 2. MERLIN images of Betelgeuse observed at 4.994 GHz in 1995 
and 1996, contours at 0.60 and 0.48 mjy beam -1 (4 cr nTls ), respectively, 
using the same angular scale as Fig.fT] The beam siz e is (85x55) mas 2 . The 
large, dashed circle shows the 310-mas, 5-GHz disc ihim et all 1993) . 

component is given by A = 7r(0 ma j# rain )/(41n2). The main compo- 
nent thus has 7b = 1 170 + 135 K. There is a point-like residual of 
0.082 mjy beam -1 (only just over 30-^) at the position of the peak, 
contributing another ~ 100 K, giving T b = 1250± 150 K within a ra- 
dius of ~5 R*. This is close to the tem perature and the size (dashed 
circle in Figs. [T] to [3} measured by iLim et af] d 1 998b by fitting a 
uniform disc to the 5-GHz VLA visibilities, showing the stability 
of these atmospheric conditions over 16 yr. 

The residuals after subtracting the central component are 
shown in the insert of Fig. [T] There is a substantial arc in the SW 
quadrant, with a maximum extent above 3cr mls of 510 mas with 
a total flux density 0.088±0.019 mjy in 0.0249 arcsec 2 . It has an 
irregular outline but is approximately 100 mas wide, at a radius be- 
tween ~ 175-275 mas from the central peak. The SW arc has an 
average T b = 150 + 40 K. 

3.2 Position and flux density 

Our first map of Betelgeuse (phase-reference solutions only) had a 
peak of 0.78±0.063 mjybeairr 1 at 05:55:10.3267 +07:24:25.525. 
The noise-based position uncertainty is 9 mas and the total astro- 
metric uncertainty (mostly du e to the 1 °5 separatio n from the phase 
reference source) is 20 mas. lHarper et al.l (J2008) used Hipparcos 
and multi-epoch VLA data to solve for position, proper motion and 
parallax, up to epoch 2004.829. The preferred solution (their num- 
ber 5) predicts a position of 05:55:10.3250 +07:24:25.536, uncer- 
tainty (5, 6) mas at our epoch, 2012.536. This is (25, 10) mas from 
our position, close to the combined uncertainties. 

The total flux density of 1.62-1.71 mjy (depending on 
whether the SW arc is enclosed) measured at 5.75 GHz, is equiv- 
alent t o (1.3-1.4)±0.2 mjy at 4.8 5 GHz, using a spectral index of 
1.32 jNewell & Hiellmind 1982). allowing for the uncertainty in 
this extrapolation. Observations at 4-8-4.9 GHz showed a decline 
from 2 to 1.2 mjy during 1981 to 2002, fluctuating within the range 



Figure 3. 1996 combined MERLIN+VLA image of Betelgeuse (contours 
at (-1, 1, 2, 4) X 0.17 mjy beam -1 ), beam size 200 mas, overlaid on the 
VLA-only image observed at 4.9 GHz, beam size 400 mas. The angular 
scale is the same as Fig. [Tl The large, dashed circle shows the 310-mas, 
5-GHz disc jLim et alj 19981) . 

0.9-1.3 mjy up to 20 04 jNewell & Hiellmind[l9^ . lskinner etal] 
ll997LlLim etal1ll998l: Harper & Brown private communication). 

3.3 Hotspots seen in 2012 

The contours in Fig.[T]show the image optimised for resolution, us- 
ing a (80 x 60) mas 2 restoring beam. Two individually-unresolved 
peaks appear, with flux densities and positions of 0.706 mjy 
beam -1 at 05:55:10.3299 +07:24:25.510 and 0.489 mjy beam" 1 
at 05:55:10.3241 +07:24:25.541, respectively (measured by fit- 
ting two Gaussian components simultaneously). cr nTls =0.0095 mjy 
beam -1 and the noise-based position uncertainties are 1 mas for the 
brighter peak and 2 mas for the fainter one. The hotspots are sepa- 
rated by 90 ± 10 mas, at a PA of (1 10 ± 10)°. This is similar to the 
orientation of the major axis of the Gaussian component fitted to 
the emission imaged at 180-mas resolution (Section [3.U and sig- 
nificantly different from the beam PA. The peaks are separated by 
> 20<r in both position and flux density. Their brightness tempera- 
tures are T b = 5400 + 600 K and T b = 3800 + 500 K, uncertainties 
dominated by the overall flux scale uncertainty. 

3.4 1995 - 1996 MERLIN and VLA results 

Figure [2] shows that, at the lower sensitivity of the 'old' MERLIN 
archive data, just a few hotspots were detected at the position of 
Betelgeuse. In 1995, 5 spots brighter than 0.60 mJy beam -1 (4o" lras ) 
were detected, the brightest at 0.78 mjy beam -1 . In 1996, 4 spots 
brighter than 0.48 mJy beam -1 (4cr rms ) were detected, the brightest 
at 0.56 mjy beam -1 . The formal position uncertainties are ~40 mas 
but may be higher due to calibration errors. The shortest baseline 
( 1 1 km) allows scales as large as 0.5 arcsec to be imaged, and the 
sensitivity of the 1995 and 1996 images corresponds approximately 
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to the highest contour in Fig. Q] so the non-detection of most of 
the stellar disc was due to lack of sensitivity rather than missing 
spacings. The peak brightness temperatures (at A = 0.061 m) were 
T b = 8000 ± 2000 K in 1995 and T b = 6000 + 2000 K in 1996. 
The uncertainties are too high to draw any firm conclusions from 
the distribution of the spots. 

The contours in Fig. [3] show the 1996 combined MER- 
LIN+VLA image, containing 1 .69±0. 17 mJy in a Gaussian compo- 
nent with (6» maJ , 9 mm ) of (289x256) mas 2 , cr g 17 mas, PA (150±30)°. 
This gives T b = 1200 + 200 K. The colour image shows the 1996 
VLA-o nly data, which w ere p art of a multi-frequency d ataset anal- 
ysed bv lLim et al.l Jl998h and lHarper & Brownl d2003l) . who mea- 
sured a flux density of 1 .77 ± 0.09 mJy. The VLA image contains a 
faint feature 0.6 arcsec due SW of the peak (over twice the distance 
to the e-MERLIN arc), at 0.095 mJy beanr 1 , just over 3cr nlis . 

4 DISCUSSION 

4.1 SW Arc 

The SW arc is at a radi us of ~250 mas, wh ere, according to the 
semi-empirical model of lHarper et al.l (2001 ). the electron temper- 
ature r e is ~900 K and the hydrogen number density «h~5x10 13 
irr 3 . This gives an optical depth ~0. 18 for the observed T b ~150 
K. Assuming that the arc is as deep as it is wide, we approxi- 
mate the region producing most of the emission by a spheroid with 
semi-major axes of 10, 10 and 25 au, giving a hydrogen mass of 
~ 1.8 x 10~ 6 M , suggesting that one such clump could be formed 
every couple of years. 

The arc is about 100 mas outside the 6-GHz p hotospheric ra- 
dius, i n a similar direction to the CN plume seen bv lKervella et al.l 
(2009) at < 6 Ri, and the faint SW extension at ~ 0.6 arcsec seen 
at 5 GHz in 1996. A c lump of CO emission was seen in 2007 by 
O' G orman et al. I J2012I) just outside the S 1 shell, 5 arcsec SW of the 
centre, suggesting a wind travel time of ~500 yr at 9 km s -1 . These 
observations suggest repeated or continuous ejections, rather than 
proper motion of a single clump. This is also aligned with t he SW 
pole of the rotational axis identified by lUitenbroek et al l dl998T) . 
The direction could b e coincidence, since observations at 7.76- 
19.5 ytim in Nov 2010 dKervelfaetai]|201lh show clumps outside 
the photosphere at 0.8-1.7 arcsec around more than a semi-circle, 
between NE and S. A comparison with the directio n of the bow 
shock seen 6-7 arcmin to the NE dDecin et al.ll2012l) is tempting, 
but structures on sub-arcmin scales are deep within the astropause 
and should share the proper motion of the star. 

4.2 Hotspots 

The 5-6 GHz hotspots, measured at 55-85 mas resolution, have 
r b formally above the peak photospheric temperature of 3600 
K lDvcketalJ[l996l) , simil ar to the chromo spheric temperature, 
which could reach 8000 K lHarper & Brownll2006l) . We used con- 
servative uncertainties, dominated by the maximum uncertainty 
in the flux scale, although the good agreement with VLA results 
for the whole-disc flux density suggests our uncertainties may be 
overestimated. One or more hotspots are significantly hotter than 
the photosphere in 2012, 1996 and 1995 and als o in the 199 2 
MERLIN observations imaged at full resolution dMorrislfeoOll) . 
Note that, at the sensitivity of old MERLIN, no emission at all 
at T b < 4500 - 6000 K (depending on epoch) would have been 
detected. 



The radio photosphere major axis is at 1 10 - 1 15° in the fully- 
calibrated 2012 data, compared with (67±7)° at 43-GHz in 1998 
dLim et alii 19981) . The orientation of the optical disc and hot spots 
also varies; the ~ 20 mas axis joining the two hotspots seen by 
Haubo is et al. I d2009l) is approxi mately NE-SW. Fo ur sets of IR 
observations over 24 months by iTufhill et alj d 19971) showed 2-3 
hots pots at a variety of posit ion angles within the inner ~ 100 mas. 

lHarper & Linskvl d200ll) showed that increasing the local T e 
2-3-fold enhances radio emission much more effectively than in- 
creasing the gas density by the same factor. There are several pos- 
sible origins for hotspots: 

(i) Levitation of gas from the inner photosphere, e.g. by convec- 
tion, although such gas would cool radiatively and by expansion. 

(ii) Exceptionally cool, small patches in the outer layers, trans- 
parent enough to expose hotter, inner layers. The separation of the 
2012 hotspots is 90 + 10 mas, similar to the 80-mas diameter of 
the hottest layer (w here the greatest chromo spheric contribution to 
emission is greatest. lllarper & Brownl2006l using 197 pc distance). 

(iii) Shock heating due to pulsations and/or convection; 

(iv) Chromospheric patches at ~5 R+ heated by acoustic and/or 
magnetic processes. 

Options (iii) an d/or (iv) seem most likely to produce radio 
hot spots >4000 K. Ilreland et all d201ll) modelled pulsation in 
Miras out to 5R+, finding average effective temperatures ^3800 
K. Large-scale up- and down-flows are seen in chromospheric 
lines, including d irectio n al reversals at up to ~3i?», measured by 
lLobel & Dupreel d200ll) . | Ohnak a et al. I dionh measured up- and 
down-draughts of CO, within ~2.5R*, at <30 km s -1 , with different 
velocities occurring simultaneously on different sides of the star. 
Such speeds are sufficient to supply shock heatin g even if it is not 
clear h ow convection could operate at »2.5i?*. Witt kowski et alj 
d201 1) made spectro-interferometric observations, around 2 fan, of 
continuum and clumpy molecular layers around AGB stars, which 
suggested pulsation- and shock-induced chaotic motion, not requir- 
ing convection. Most radio emission is non-chromospheric (Sec- 
tion [T|l i _birtJhe_jiotepots have Tb <8000 K, in the chromospheric 
range. Harp er & Brownl d2006l) suggest that chromospheric gas is 
confined within magnetic fields; dissipation of magnetic energy 
may continue to heat the gas even at high elevations. 



5 FURTHER HIGH-RESOLUTION STUDIES 

We have shown that both cool, extended radio-continuum emis- 
sion and hot starspots on Betelgeuse can be resolved at 5.5-6 GHz. 
Multi-epoch, multi-frequency follow-up is needed to establish the 
emission mechanisms and role in mass loss. In full operations, e- 
MERLIN will reach 0.002 mJy beam -1 in two full tracks at 5-7 
GHz. Proper motions of 5 km s could be measured in 4 months 
for a spot 0.3 mJy brighter than its surroundings, revealing whether 
the hotspots rotate as predicted for the chromosphere. The vari- 
ability ti mescale will constrain t he underlying mechanism. Mod- 
elling bv lHarper & Linskvl J200ll) shows that CII recombination to 
CI takes ~ 2 months at ;jh~4x10 14 irr 3 , around 5 R*, or longer at 
higher radius, and chemical changes have similar timescales, but 
mechanisms involving bulk up-/down-draughts alone would take 
about a year to produce comparable changes. 

Matching resolution at 50 mas can be achieved by the EVN, 
e-MERLIN and the Karl G. Jansky VLA (alone or in combination 
as required) at 1.4-50 GHz, and eventually by ALMA at higher 
frequencies. Finer resolution (< 20 mas) is possible at 22 GHz and 
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some higher frequencies, suited to the smaller observed size of the 
radio photosphere. This will allow us to resolve the spectral indices 
of the extended emission and the ho tspots and test models for lo- 
calised heating dHarper et all d200ll) and references therein). For 
example, a decrease of T s with decreasing frequency above a hot 
spot would suggest a chr omospheric, rath er than convective origin. 

The average radius dLim et alii 19981) is approximately related 
to the observing frequency in GHz by R(v) ~ 55v~ 05 au. At a speed 
of 10 km s" 1 (within th e range shown by chromospheric lines, 
iGilliland & Dupredl 19961) . it would take 0.5 yr for a disturbance to 
propagate 1 au from surfaces optically thick at 25 GHz to 21 GHz, 
and 1 yr from the 6 GHz to the 5 GHz surface. Monitoring at de- 
creasing frequencies, i.e. increasing radii, at suitable intervals, will 
track the evolution of disturbances. This will show whether we are 
witnessing phenomena propagating outwards smoothly, or local- 
ized changes behaving differently, helping to distinguish between 
convection and pulsation. We will resolve CO, other molecules and 
dust using ALMA, thus measuring the composition and mass of the 
wind. Clumps will be resolved on scales of a few au, within a few 
tens au from the star, where chemi stry is still active, ov erlapping 
with the scales resolved in the MIR dKervella et alj[20llh . 



6 SUMMARY 

e-MERLIN has produced the highest-resolution images of Betel- 
geuse at cm wavelengths, in July 2012. The 5.75 GHz photosphere 
is fitted by a 2-D Gaussian component, axes (0.235x0.218) arcsec 2 
(~ 5i?*), flux density 1.62 mjy, with an average brightness temper- 
ature ~ 1250 K. Extrapolating to 4.9 GHz shows that the flux den- 
sity decline since 1981 may have levelled off; t he large-scale siz e 
and brightness is consistent with the results from lLim et al. I dl998l) . 
The peak position is close to the positio n predicted from the proper 
motion analysis bv lHarper et al.l 112008), An arc is seen 0.2-0.3 arc- 
sec to the SW, brightness temperature ~ 150 K, hydrogen mass 
~ 1.8 x 10~ 6 M . It is in the same direction as IR plumes seen 
on sm aller and larger scales dKervella et al.l [20091 : iKervella et al.l 
l201ll) . a faint 4.9-G Hz feature at 0.6 arcsec seen in 1996 and CO 
emission imag ed bv lO'Gorman et al.l d2012l) . This suggests persis- 
tent ejection in this direction for several centuries. 

e-MERLIN resolved two distinct peaks, 90 mas apart within 
the larger stellar disc, with brightness temperatures 4000-6000 K. 
These are likely to be hot spots at ~5 R+, although it is possible 
that cooler, low-opacity inhomogeneities in higher layers reveal 
the hottest chromospheric layers. Similar, isolated peaks were also 
detected by MERLIN observations in 1992, 1995 and 1996. Such 
high-sensitivity, high-resolution radio continuum observations can 
resolve emission simultaneously from material at chromospheric 
and photospheric temperatures, and have the potential to track 
proper motions. Future multi-frequency imaging with e-MERLIN, 
the VLA and ALMA will investigate whether the signatures of 
chromospheric emission, convection or pulsation shocks are seen 
in the radio photosphere from 2-10 i?*, helping to elucidate the 
mechanisms which eject material from the star into its envelope. 
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APPENDIX A: DATA PROCESSING SUMMARY 
Al Calibration 

This is a summary of the steps for data calibration. AIPS dGreisenl 
1 1994b was used for data processing as the required capabilities are 
not yet fully implemented in CASA. Data were stored in multi- 
source format, the time-dependent solution tables being interpo- 
lated into calibration tables using CLCAL. Separate solutions were 
derived for each spectral window and hand of circular polarization. 

(i) Load data (FITLD), ensure that they are in time-baseline order 
(MSSORT) and inspect/edit. Editing (SPFLG, UVFLG and IBLED) 
was performed at a number of stages as required. 

(ii) Derive delay and rate corrections per scan for calibration 
sources 0551+0829 and OQ208 (FRING). 

(iii) Derive time-dependent phase, and then amplitude and phase 
solutions for calibration sources (CALIB). 

(iv) Use the solutions for the three most similar anten- 
nas (Darnhall, Pickmere, Knockin) for a selected time 
interval, to establish the flux density of 0551+0829 by 
comparison with OQ208 (GETJY), refining the solu- 
tions by requiring a straight-line spectral index (SOUSPC). 

Centre freq. (GHz) 5.4954 5.6874 5.8154 5.9434 
OQ208 (Jy) 2.48 2.49 2.50 2.51 

0551+0829 (Jy) 0.183 0.190 0.196 0.203 

The flux density of OQ208 was derived from 3C286 (R. Beswick, 

private communication). 

(v) Derive channel-by-channel amplitude and phase bandpass 
solutions for OQ208, averaging each half-hour scan (BPASS). 

(vi) Apply the bandpass solutions using SPLAT and repeat steps 
(ii), (iii) and (v), using the established flux density of 0551+0829 
as its model. 

(vii) Split out Betelgeuse, applying all calibration. 

(viii) Make a preliminary, naturally-weighted image (IMAGR). 
The star was clearly resolved, showing a bright central ridge and 
a halo a few hundred mas in size, strongest to the SW. AIPS fit- 
ted beam size of (0.195x0.074) arcsec, PA 143°. The major axis is 
more than double the true resolution (even at declination 7°), be- 
cause the great sensitivity of the Lovell compared to the other tele- 
scopes creates a bias towards shorter baselines. All imaging used 
careful boxing and a small number of iterations (50-100). 

(ix) Use the clean components as a model for phase-only self- 
calibration with a 7-min (per-scan) solution interval. These were 
applied and a new image of Betelgeuse was made. We also plot- 
ted the visibility amplitudes and noticed some residual, slowly- 
changing amplitude errors (lower amplitudes on shorter baselines, 
differences between hands of polarization). We carried out a final 
amplitude self-calibration with a 1 hr solution interval, which did 
not affect the source structure but improved the signal to noise ratio. 

A2 Imaging 

All images were made in total intensity, combining all spectral win- 
dows using a spectral index of 1.32 (using a flat spectral index 
made a very small difference, much less than the 10% total flux 
scale uncertainty). There is a trade-off between sensitivity to ex- 
tended structure and resolution, in interferometry data. We chose 
two weighting schemes, based on the inherently different sensitiv- 
ities of the e-MERLIN telescopes and on empirically established 
values to optimise either sensitivity or resolution without increasing 
the noise more than 10%. The Defford telescope is least sensitive 
and was given a weight of 0.1. The 75-m Lovell telescope (provid- 



ing many short baselines) and the 32-m Cambridge telescope (pro- 
viding the longest baselines) are the most sensitive at 5-6 GHz. We 
maximised sensitivity to extended structure by giving Lovell and 
Cambridge data weights of 10 and 2.4, respectively. We maximised 
resolution by giving both the Lovell and Cambridge data a weight 
of 4. All other (25-m) telescopes were given a weight of 1. The 
weighting was further modified in the AIPS task IMAGR by apply- 
ing suitable interpolation into the missing spacings; a uvtaper of 
3500 kA with natural weighting for optimum sensitivity and partial 
uniform weighting ROBUST 0.75 for optimum resolution. 

A3 Astrometry and photometry 

We observed Betelgeuse at a pointing position of J2000 
05:55:10.33 +07:24:25.6. The phase reference ICRF2 position of 
05:51:11.2293 +08:29:11.221 is accurate to < 0.001 arcsec and it 
is only l.°5 from Betelgeuse. Inspection of the OQ208 raw phase 
corrections suggests this introduces an uncertainty of one turn of 
phase in 30 min of time on the longest baselines, corresponding 
to one beam in 7°5 of arc or 0.016 arcsec in l.°5; the uncertain- 
ties in telescope positions introduce negligible errors. We measured 
the position of Betelgeuse for astrometric purposes by fitting a 2D 
Gaussian component to the naturally-weighted image before self- 
calibration. 

The flux density scale is derived from 3C286; we do not yet 
have a full model of this source at e-MERLIN resolution and sen- 
sitivity so we conservatively estimate a 10% flux scale uncertainty. 
0551+0829 has a spectral index of 1.3. We note that MERLIN ob- 
servations in 1995 found a flux density of 0.128 Jy for 0551+0829 
at 4.994 GHz, 20% lower than the flux density extrapolated from 
current observations. This is within the plausible variability of a 
compact QSO in 17 yr. 
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